clc;
syms x y z;
f(x,y,z)=sin(x^2*y)*exp(-x^2*y-z^2);
df=diff(f,x,x,y,z);%对x求2次偏导
F=simplify(df)
pretty(F)
latex(F)